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^ ; Abstract 
m ■ 

o : 

O ; 

' We study the fluctuation properties of AT = electromagnetic transition 

' 

' intensities and electromagnetic moments in ^ ~ 60 nuclei within the frame- 

^3 : 

^ ' work of the interacting shell model, using a realistic effective interaction for 

> : 

' pf-shell nuclei with a Ni core. The distributions of the transition intensities 

X; 

■ and of the electromagnetic moments are well described by the Gaussian or- 
thogonal ensemble of random matrices. In particular, the transition intensity 
distributions follow a Porter-Thomas distribution. When diagonal matrix ele- 
ments (i.e., moments) are included in the analysis of transition intensities, we 
flnd that the distributions remain Porter-Thomas except for the isoscalar Ml. 
The latter deviation is explained in terms of the structure of the isoscalar Ml 
operator. 
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1 



Typeset using REVT^ 



Random matrix theory (RMT) [jT| was originally introduced to explain the statistical 
fluctuations of neutron resonances in compound nuclei 0. The theory assumes that the 
nuclear Hamiltonian belongs to an ensemble of random matrices that are consistent with the 
fundamental symmetries of the system. In particular, since the nuclear interaction preserves 
time-reversal symmetry, the relevant ensemble is the Gaussian orthogonal ensemble (GOE). 
The use of RMT in the compound nucleus was justified by the complexity of the nuclear 
system. Bohigas et al |^ conjectured that RMT describes the statistical fluctuations of 
a quantum system whose associated classical dynamics is chaotic. RMT has become a 
standard tool for analyzing the universal statistical fluctuations in chaotic systems 0-0. 

The chaotic nature of the single-particle dynamics in the nucleus can be studied in 
the mean-field approximation. The interplay between shell structure and fluctuations in the 
single-particle spectrum has been understood in terms of the classical dynamics of the nucle- 
ons in the corresponding deformed mean-field potential . However, the residual nuclear 
interaction mixes different mean-field configurations and affects the statistical fluctuations 
of the many-particle spectrum and wave functions. These fluctuations can be studied using 
various nuclear structure models. The statistics of the low-lying collective part of the nu- 
clear spectrum have been studied in the framework of the interacting boson model [^jjl^], in 
which the nuclear fermionic space is mapped onto a much smaller space of bosonic degrees 
of freedom. Because of the relatively small number of degrees of freedom in this model, it 
was also possible to relate the statistics to the underlying mean-field collective dynamics. At 
higher excitations, additional degrees of freedom (such as broken pairs) become important 



111 , and the effects of interactions on the statistics must be studied in larger model spaces. 
The interacting shell model offers an attractive framework for such studies, in which realistic 
effective interactions are available and where the basis states are labeled by exact quantum 
numbers of angular momentum, isospin and parity [0 . 

RMT makes definite predictions for the statistical fluctuations of both the eigenfunctions 
and the spectrum. The electromagnetic transition intensities in a nucleus are observables 
that are sensitive to the wave functions, and the study of their statistical distributions should 



complement @JT^ the more common spectral analysis. Using the interacting shell model, 
B{M1) and B{E2) transitions were recently analyzed in ^^Na |jl3[ and found to follow the 



Porter-Thomas distribution WM, in agreement with RMT and with the previous finding of 



a Gaussian distribution for the eigenvector components |T^- . 

Most studies of statistical fluctuations in the shell model have been restricted to lighter 
nuclei {A ^ 40) where complete Ohu calculations are feasible (e.g., sd-shell nuclei). Here 
we study the fluctuation properties of AT = electromagnetic transition intensities in 
nuclei with A ~ 60. We also study the statistics of the diagonal matrix elements, which 
describe the corresponding electric or magnetic moments. We find that the statistics of both 
the transition intensities and the electromagnetic moments are well described by RMT. In 
RMT, off-diagonal matrix elements have a Gaussian distribution with zero mean, leading 
to a Porter-Thomas distribution of the transition intensities. Diagonal elements also have 
a Gaussian distribution but with non-zero mean. Thus the squares of the diagonal matrix 
elements should also follow a Porter-Thomas distribution but only after the mean value is 
subtracted from the diagonal elements. Experimentally it is difficult to measure the moments 
of excited states, while transition intensities are more readily accessible. 

When the good quantum numbers (isospin, spin and parity) of the initial and final states 
are the same, we can include the diagonal matrix elements in the analysis of the transition 
intensity distribution. The number of diagonal elements is much smaller than the number 
of off-diagonal elements, and generically we do not expect the diagonal elements to modify 
the Porter-Thomas statistics of the transition intensities. This is confirmed in our shell 
model calculations with the exception of the Ml operator in self-conjugate nuclei (T^ = 0), 
where a significant deviation from a Porter-Thomas distribution is found.[| This deviation 
is explained in terms of the special structure of the isoscalar Ml operator, which leads to 
diagonal matrix elements that are much larger than its off-diagonal elements. 



We find a similar deviation in sd-shell nuclei. 
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Full p/-sliell calculations are currently not feasible for A ~ 60 nuclei. We therefore 
perform p/-sliell calculations with ^^Ni as a core, i.e., we assume fully occupied /7/2 orbits 
and consider all possible many-nucleon configurations defined by the O/5/2, IP3/2 and lpi/2 
orbitals. The effective interaction is chosen to be the isospin-conserving F5P interaction ||19|| . 
This interaction is successful in describing the mass range A ~ 57 — 68. The calculations 



were performed using the shell model program OXBASH pO . 

To test the validity of RMT in the above model space, we first study the spectral fluc- 
tuations for states with good spin and isospin (all states in the p/-shell model space have 
the same parity). Fig. |l|(a) shows the nearest-neighbors level spacing distribution P{s) 
for the unfolded T = levels. To improve the statistics we combine the spacings of the 
J = 2 states and the J = 4 states. The calculated distribution (histograms) is in agreement 
with the Wigner-Dyson distribution (solid line), and level repulsion is clearly observed at 
small spacings. The Poisson distribution, which corresponds to a random sequence of levels 
(and characterizes regular systems), is shown by the dashed line for comparison. Another 
measure of fluctuations is the spectral rigidity described by Dyson's A3 statistic The 
bottom panels of Fig. Q show A3(L) versus L for J = 2 states with T = [panel (b)] and 
T = 1 [panel (c)]. The results agree well with the GOE limit (solid lines). 

We now consider the electromagnetic transition intensities. Denoting by B(pL;i — > /) 
the reduced transition probability from an initial state \i) to a final state |/), with u indicat- 
ing the electric (E) or magnetic (M) character of the transition, and 2^ the multipolarity, 



we have [21 



B{uL: Jjlilz — > Jf-L fiz) = ; , , ^ ^ . (1) 

^ ' ^ * ^ ^ ^ (2Ji + l)(2ri + l) ^ ^ 

Here MisipL) and MmipL) are the triply reduced matrix elements for the isoscalar and 
isovector components of the transition operator, respectively. Note that these matrix ele- 
ments depend on Jj, Tj and J/, Tj but not on T^. For AT = transitions (Tj = Tf = T) the 
isospin Clebsch-Gordan coefficient in Eq. (||) is simply given by 



{TTzW\TTz) = Tz/JT{T + 1) . (2) 



It follows that the isovector component in Eq. is absent for self-conjugate (i.e. Tj = 0) 
nuclei. Thus the statistics of the isoscalar component of an electromagnetic transition op- 
erator can be inferred directly from the study of AT = transitions in = nuclei. 
Consequently, we are able to test the sensitivity of the statistics to the isovector and the 
isoscalar contributions. Below we present results for E2 and Ml transitions. The E2 tran- 
sitions were calculated using standard effective charges of = 1.5 and e„ = 0.5. We note, 
however, that the isoscalar and isovector components of E2 are (up to a proportionality 
constant) independent of the effective charges, and their corresponding statistics are thus 
also independent of the particular choice of effective charges. 

To study the fluctuation properties of the transition rates, it is necessary to divide out 
any secular variation of the average strength function versus the initial and final energies. 
We do this by applying the method of Refs. and |10[ . We calculate an average transition 



strength at an initial energy E and final energy E' from 

Ei f B(pL; I /)e-(^-^0V27^e-(s'-i?/)V27^ 
{B{ujL- E, E')) = ^^je-(i^-i^0V27^e-(^'-^/)V27^ ' 

where 7 is a parameter chosen as described below. For fixed values of the initial (</«, T) and 

final {Jf,T), we calculate from Eq. ([I|) the intensities B{ujL;i /). All transitions of a 

given operator (e.g.. Ml or E2) between the initial and final states of the given spin and 

isospin classes have been included in the statistics. We remark that the energy levels used 



in (^) are the unfolded energy levels characterized by a constant mean spacing. The 
value of 7 in (|^) has been chosen to be large enough to minimize effects arising from the 
local fluctuations in the transition strength but not so large as to wash away the secular 
energy variation of the average intensity. In the present calculations we use 7 ~ 2.5. We 

renormalize the actual intensities by dividing out their smooth part 

^ BjZJL; J{TT, ^ J^TT,) 

{BipL-E,E')) ' ^' 

and construct their distribution using bins that are equally spaced in log^g y- The choice 
of lo 

gio y the variable allows us to display the distribution of the weak transitions over 



several orders of magnitude p4 



In RMT we expect a Porter-Thomas distribution for P{y), i.e., a distribution distri- 
bution in z/ = 1 degrees of freedom [^]. A distribution in u degrees of freedom is given 
by 

= {u/2 < y >rl^y^l^-^e-^yl^<y>lV{yl2) . (5) 

Indeed, consider a transition operator T and its matrix elements (/|T|i) between an initial 
state \i) and a final state |/). Using a complete (and fixed) set of states |m), we have 

(/|rK) = ^r„,„^*(m)7A,H, (6) 

m,n 

where '^i[n) = {n\i) represents the eigenfunction \i) in the fixed chosen basis and Tm,n = 
{m\T\n). In RMT, V^j(n) are Gaussian random variables, and by the central limit theorem 
it follows from that {f\T\i) is also Gaussian. In RMT, two different eigenfunctions 



are uncorrelated il)*f{m)il)i{n) = 0, and consequently {f\T\i) = 0. Similarly, using Wick's 
theorem for Gaussian variables, the variance is given by 



iifin)) = E T:^,nTm',n'M^)rf{m') Mn)^:in') = ^Tr (P.T^PfT) , (7) 

7X1,71,171' ,n' * / 

where Pi is the projection operator on the A^j-dimensional subspace of eigenstates with quan- 
tum numbers Ji\Ti {Pf is similarly the projection operator on the corresponding subspace 
of final states). The transition intensity is proportional to the square of the matrix element, 
and thus has a Porter-Thomas distribution (i.e., Eq. (^ with u = 1). 

We examined E2 and Ml 2+,T = 1 — 2+,T = 1 transitions in A = 60 nuclei. The 
T = 1 states form isobaric multiplets with = 0, ±1 (i.e., ^°Co, ^°Zn and ^°Cu). We 
studied the statistics of the transition intensities in both = and = I nuclei. Because 
of the vanishing of the isovector Clebsch-Gordan coefficient in Eq. (|1|), the transitions in 
the self-conjugate nucleus ^°Zn (T^ = 0) are purely isoscalar. For each transition opera- 
tor we sampled 66^ — 66 = 4290 transition matrix elements. The calculated distributions 
(histograms) of the B{E2) (left panels) and B{M1) (right panels) 2+, T = 1 2+, T = 1 
transitions are shown in Fig. ^. The top and bottom panels correspond to Tj; = 1 and = 

6 



nuclei, respectively. All distributions are well described by the Porter- Thomas distribution. 
Similar results were found by Adams et al |T3] for the N = Z (T^ = 0) nucleus ^^Na in the 



sci-shell. 

When the quantum numbers of the initial and final states are identical, we can also 
examine the statistics of the diagonal matrix elements {i\T\i). Using with i = f and the 
central limit theorem we still find that {i\T\i) is a Gaussian variable but with a non-zero 



average 
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{t\T\t) = J2 %n,ni^tMUn) = ^Tr (P.T) . (8) 

m,n i 

Using Wick's theorem and the fact that the GOE wave functions ipi{m) are real, the variance 
of a diagonal element is 



For a real and hermitean operator T, we find 

<y'{m^) = (P.TP.T) , (10) 

which is a factor of two larger than the variance of an off-diagonal element (see Eq. (|^ for 
the case when the initial and final states span the same subspace of states, i.e.. Pi = Pf). 

Since the mean value of {i\T\i) is non-zero, its square does not follow a Porter- Thomas 
distribution. However, once the mean value (H) has been subtracted from the diagonal 
matrix elements, their squares z are predicted to have a Porter-Thomas distribution. The 
distributions of the squares of these shifted reduced diagonal matrix elements are shown in 
Fig. ^ for the Ml magnetic moments of the J = 2+,T = 1 states in A = 60 nuclei. The 
values of z are renormalized by dividing out the secular variation with energy of the square 
of the reduced diagonal matrix element, defined by an equation similar to (^ but with a 
single Gaussian (since i = f). Considering the small number of data points used to compute 
the distribution P(log]^Q z) (there are only 66 diagonal matrix elements), the agreement with 
a Porter-Thomas distribution (solid lines) is reasonable. 



It is of interest to see the effect that the inclusion of diagonal matrix elements has on the 
distribution of transition intensities. Since the number of diagonal elements is much smaller 
than the number of transition intensities, we do not expect them to modify the Porter- 
Thomas statistics of the transition intensities (even though the squares of the unshifted 
diagonal reduced elements do not follow a Porter-Thomas distribution). Fig. ^ shows the 
same distributions as in Fig. |^ but when the squares of the reduced diagonal elements of 
the corresponding transition operator are included in the analysis. The B{E2) distributions 
as well as the B{M1) distribution in the Tj = 1 nucleus are all well described by a Porter- 
Thomas distribution, but the B{M1) distribution in the self-conjugate nucleus {Tz = 0) 
deviates significantly from the Porter-Thomas distribution (right bottom panel of Fig. 
The Ml transitions in self-conjugate nuclei (T^ = 0) have the property that they are purely 
isoscalar. In Tj = 1 nuclei both isoscalar and isovector components contribute to the Ml 
transitions, but the isoscalar Ml matrix elements are much weaker than the isovector Ml 
|2T| , and the latter dominate the distributions. To check that our conclusion does not depend 
on the size of the model space, we performed similar calculations in a model space with 6 
active particles, i.e., for A = 62 nuclei. Again, we observed that the distributions of the Ml 
transitions in the self-conjugate nucleus = (i.e., ^^Ga) deviate from a Porter- Thomas 
distribution. 

The deviation from a Porter- Thomas distribution for the isoscalar Ml can be traced to 
the diagonal matrix elements. Although their number is much smaller than the number of 
off-diagonal elements, they are in magnitude much larger than the off-diagonal elements. 
Indeed, the isoscalar Ml has a particularly simple structure and can be written (in units of 
nuclear magneton fi^) as 

^L + ^^S, (11) 

where L and S are the total orbital angular momentum and spin, respectively, and g^ .^, 
Qp^ are the proton and neutron orbital and spin g factors. Since the transition matrix 
elements of the total angular momentum J = L + S between different states vanish, the 
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orbital and spin contributions to the isoscalar Ml matrix elements always have opposite 
signs. Consequently, the off-diagonal elements are small. However the interference between 
the orbital and spin contributions to the diagonal matrix elements can be constructive, 
leading to larger diagonal matrix elements. These diagonal matrix elements are related to 
the magnetic moments, which in = nuclei are given by J/2 + 0.38 (S^) PT|. The mean 
value (H) of the isoscalar magnetic moment is large (relative to a typical off-diagonal element) 
and leads to the observed deviation from a Porter-Thomas distribution. In the case of the 
isovector Ml matrix elements, the relative phase of the orbital and spin contributions is 
more random and the corresponding statistics are Porter-Thomas even when diagonal 
elements are included. 

In conclusion, we have studied AT = B{E2) and B{M1) transition strength distribu- 
tions in p/-shell nuclei with A ~ 60 and found that they agree well with Porter-Thomas 
statistics. We have also studied the statistics of diagonal elements (e.g.. Ml magnetic mo- 
ments) in these nuclei. The distribution of the squares of these moments is found to be in fair 
agreement with a Porter-Thomas distribution, but only when the mean value is subtracted 
from the moments. In general, the inclusion of off-diagonal elements in the statistics of tran- 
sition intensities does not affect the observed Porter-Thomas statistics. An exception is the 
isoscalar Ml operator whose diagonal matrix elements are much larger than its off-diagonal 
elements. 

This work was supported in part by the Department of Energy grant No. DE-FG-0291- 
ER-40608 and by the Russian Foundation for Basic Research, project No. 00-02-17194. 
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FIG. 1. Spectral statistics in ^4 = 60 nuclei, (a) Nearest-neighbor level spacing distribution 
for r = states with J = 2,4 (histograms). The solid line is the Wigner-Dyson distribution 
and the dashed line is the Poisson distribution. Bottom: A3 statistic for J = 2 levels with (b) 
T = and (c) T = 1. Also shown are the GOE (solid lines) and Poisson (dashed lines) limits. 
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FIG. 2. The B{E2) (left panels) and B{M\) (right panels) intensity distributions (histograms) 
for the 2^^ ,T = 1 ^ 2'^,T = 1 transitions in ^ = 60 nuclei: (a,b) = 1; (c,d) = 0. The solid 
lines describe the Porter-Thomas distribution (Eq. (^) with = 1). 
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FIG. 3. Statistics of the Ml magnetic moments of the J = 2+,T = 1 states in yl = 60 nuclei: 
(a) Tz = 1; (b) Tz = 0. The histograms show the distribution P(logxo-z) where z are the squares 
of the Ml diagonal reduced matrix elements whose mean value (^) has been subtracted. The solid 
lines describe Porter-Thomas distributions. 
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FIG. 4. B{E2) (left panels) and B{M1) (right panels) intensity distributions (histograms) for 
A = 60 nuclei when diagonal matrix elements are included in the analysis. Shown are the same cases 
as in Fig. |2[ Notice that the B{M1) distribution deviates from the Porter-Thomas distribution in 
the self-conjugate = nucleus. 
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